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I. INTRODUCTION 



The Einstein gravitational field equations with positive cosmological constant [l( have 
a highly symmetric matter- free solution, de Sitter spacetime [1, 0] . Nearly a century after 
the discovery of this mathematical solution, de Sitter spacetime occupies a central place in 
modern theoretical physics and observational cosmology (see, for example, the reviews ^B])- 
It is, then, all the more interesting if something new can be said about de Sitter spacetime, 
even if the context is nonstandard. 

In recent work on the cosmological constant problem, we noted parenthetically (Foot- 
note 1 in Appendix A of Ref. that, for the simple model considered, de Sitter spacetime 
corresponded to an unstable critical point. The simple model considered 0,0] had a classical 
vector field V a (x) with a "wrong-sign" kinetic term [giving energy density p vec < for the 
cosmological solution], which we suspected to be responsible for the de Sitter instability. It 
will, however, be shown in the present article that the de Sitter instability is also present in 
the model with a "correct-sign" kinetic term [giving p vec > for the cosmological solution]. 

The particular type of vector-field theory considered (Sec. |TT]) is, most likely, pathological, 
having instabilities at the classical level and ghosts at the quantum level. Still, the vector 
field interacts only gravitationally with the other matter fields. As such, this classical 
vector field may be used to describe certain nonstandard gravitational effects in the long- 
wavelength (low-energy) limit. Two examples of such effects are discussed in the present 
article, namely, a particular type of instability of the de Sitter equilibrium solution (Sec. IHIj) 
and the corresponding final singularity (Sec. lIVp . The estimated de Sitter decay time and the 
violation of certain energy conditions by the relevant vector-field configurations are discussed 
in Sec. |V] In that last section, it is also explained how this type of classical vector field can 
perhaps play a role in the macroscopic description of a fundamental quantum instability of 
de Sitter spacetime [9-11 . 



II. THEORY 

Consider general relativity with a positive cosmological constant A and a single classical 
vector field V a (x). The specific gravitational model p-|8( used here has the following action 
(c = h=l): 



0] , (la) 



S[g, V, 0] = - j d 4 x y/-det(g) Q (Planck) 2 R\g\ + e(Q 1 [g, V})+A + £ M [g, 
e(Qi[g, V}) = -(Q 1 [g, V}) 2 = -V^V"* , (lb) 

^Planck = (87rG)- 1/2 , G > , A > , (lc) 

where a generic massless matter field <f>(x) has been added with a standard Lagrange density 
Cm(x). The action flTaJ is really classical, but, for convenience, we use quantum terminology 
such as -Epianck- I n principle, it is also possible to add a mass term for the vector field, but 
we refrain from doing so for the moment and the theory maintains the shift invariance of 
the vector field. 
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Notice that, unlike the case of a gauge field with a Maxwell action-density term, the time 
derivative of the Vq component enters the action-density term (fibl) . It is, of course, known 
that, in Minkowski spacetime (A = 0), gauge invariance is re quir ed for the Poincare invari- 



ance, locality, and stability of the massless- vector-field theory [12]. However, as explained in 
Sec. HI our interest in the classical massless vector field from ([1]) is only as an effective way 
to describe possible nonstandard gravitational effects related to the cosmological constant 
A > 0. Our focus will be on stability issues in a cosmological context. 

Let us restrict our attention to the spatially flat (k = 0) Robertson- Walker metric [3] with 
a perfect-fluid standard-matter component and an isotropic vector field (vanishing spatial 
components in appropriate coordinates). The dimensionless ordinary differential equations 
(ODEs) are then 

3h 2 = i + {j 2 + 3h 2 v 2 + X^tm , (2a) 
2h = 2hv 2 + Ahvv -2v 2 - X' 1 {l + w M )r M , (2b) 
d + 3hv-3h 2 v = 0, (2c) 
r M + 3(l + w M )hr M = 0, (2d) 

where a numerical factor \f\ = V / A/(-Epi an ck) 2 > has been absorbed into the definitions 
of the dimensionless inverse Hubble parameter h~ l and the dimensionless cosmic time r 
(the overdot stands for differentiation with respect to this r). The dimensionless variable v 
corresponds to the vector-field time-component Vo and the dimensionless variable tm > 
corresponds to the standard-matter energy density pu > with constant equation-of-state 
parameter % > 0. See Appendix A of Ref. j(| for further details. 

Using Eq. (A6) from Ref. jsj, the corresponding dimensionless vector- field energy density 
r vec and pressure p vec are found to be given by 

r V ec = e(q 1 )-q 1 ^- = +(q 1 ) 2 = \(v 2 + 3h 2 v 2 )>0, (3a) 
dqi 

Pvec = — r vcc — 2\{hv 2 + 2hvv — v 2 ) , (3b) 

parts of which, divided by A, can be seen to appear on the right-hand sides of fl2"a|) and 
( 12bj) . Remark that, in a Minkowski background with H = A = 0, the vector-field fluid (J3]) 
is not unusual, it corresponds to a matter component with an ultrahard equation of state, 
p vec = -Pvec = (dVo/dt) 2 > 0. What is unusual is how the vector-field pressure f )3b|) behaves 
in a nonflat spacetime background, possibly having P vcc < —p vcc < 0. 

III. UNSTABLE EQUILIBRIUM 

The ODEs (F2]) have an asymptotic equilibrium solution (critical point) corresponding to 
de Sitter spacetime: 

l/h(r) = V3, (4a) 
v(t) = 0, (4b) 
r M {r) = 0. (4c) 
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FIG. 1: Numerical solutions of ODEs ((2b]) and ([2c]) with dimensionless cosmological constant A > 
and vanishing standard-matter component, tm(j) = 0. The left panel shows the dimensionless 
vector-field component v(t). The middle panel shows (3/i 2 ) -1 v 2 /(l — v 2 ). The right panel shows 
^vec/(3A/i 2 ) as the ascending curves and 3/4 + (1/4) (p ve c/r vec ) as the descending curves. The 
vector-field boundary conditions are v(0) = 1/130 and v(0) = ±1/130, where the solid and dashed 
curves correspond to the plus and minus sign, respectively. The values for h(0) follow from (I2ap 
[both curves have, in fact, the same value of h(0)]. The Hubble parameter h{r) and the vector-field 
energy density r vec (r) diverge at r w 9.635 for the v(0) > boundary condition (ascending solid 
curve of the right panel) and at r ~ 11.794 for the v(0) < boundary condition (ascending dashed 
curve of the right panel) . 



There are no asymptotic solutions with h ~ 0, which would approach Minkowski space- 
time (g| [s|] . The explanation is that, with the minus sign chosen in ( ITb|) . the cosmological 
constant cannot be canceled by the vector-field contribution, as the right-hand side of ( l2aj) 
makes clear {\~ 1 tm is non-negative). 

The equilibrium solution (j4]) is, however, unstable. We will show this numerically, but it 
can also be proven mathematically by following the discussion in Appendix A of Ref. jf| (in 



fact, the linearized analysis suffices, according to Theorem 3.2 of Ref. |13|). 

Numerical solutions have been obtained with vanishing and nonvanishing standard-matter 
components. For the case of a nonvanishing standard-matter component, it is found that an 
asymptotic de Sitter spacetime is approached if the vector field is strictly equal to zero, but 
not if the vector field is nonzero. As the conclusion is essentially the same for the case of a 



vanishing standard-matter component, we focus on the tm = case 14j. Instead of reaching 



an asymptotic de Sitter spacetime, the model universe of Fig. [T]is seen to terminate after a 
finite time interval, having a diverging Ricci scalar R oc —6 (h + 2 h 2 ) at r ~ 10. The same 
type of behavior as shown in Fig. [1] is obtained for boundary conditions taking values in the 
finite intervals v(0) G (0, 1/100] and v(0) G [-1/100, +1/100]. 

Hence, we have established numerically the classical instability of de Sitter spacetime 
01]) under vector-field perturbations Sv(t) which break the original de Sitter symmetry. We 
have, in addition, explicit analytic results for the linear perturbations but will not present 
them here, as the numerical results suffice to demonstrate the instability. 



IV. FINAL SINGULARITY 

The numerical results of the previous section suggest that the model universe of Fig. [T]runs 
into a final singularity (also known as a big-rip-type future singularity or, more generally, 
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as an exotic future singularity, see Refs. 15M22| and references therein). Some analytic 
results have been obtained for the vector-field theory (TjQ) with vanishing standard-matter 
component, tm(t) = 0. The pure vector-field theory considered here is of interest, because 
it has a strictly non- negative energy density ( |3a|) . different from the scalar-field theory with 



a negative quartic coupling constant as discussed in Ref. [19 
The ODEs (TSl) now reduce to 



(hy 



v + 3hv — 3h 2 v 



2 h vv — v £ 







hr 2 + hr 2 ij 2 
3(l-v 2 ) 



(5a) 
(5b) 

(5c) 



where h, h, and 1 — v 2 are assumed to be nonzero. Next, make a change of variable s = ln(a) 
for cosmic scale factor a(r), with a(0) = 1 and Hubble parameter h = a/ a = dsjdr. The 
following ODEs for v(s) and h(s) are found: 

2v-v' 



v + 



(v'f + 3 (v' - v) 



(h>) 



1-v 2 
_i i2v-v' 







1-v 2 



v'h 



h 



-2 



-1 



,A2 



3(1 



(6a) 
(6b) 
(6c) 



where the prime stands for differentiation with respect to s. There are only two arbitrary 
constants of integration as the last two ODEs in Jp are first-order (the first ODE is consistent 
with the last two ODEs; cf. Sec. Ill A of Ref. fl). 

It is easy to check that a particular combination of trigonometric functions solves the 
nonlinear ODE (|6a|) . With an arbitrary real amplitude A G [—1, +1] and an arbitrary 
relative sign entering the solution for v(s) and with a nonzero real amplitude B in h(s), the 
following solutions of the ODEs f l6a|) and (I6bj) are obtained: 



v(s) = A sin (VS s) ± Vl-A 2 cos (VZ s) , (7a) 

l/h{s) = B [1 + (2A 2 - 1) cos^v 7 ^ s) t 2A Vl-A 2 sin(2V3 s)] exp(3 s) . (7b) 

A consistent solution of the differential system ([6]) with v and h given by (17a|) and (17bj) 
requires a solution of the constraint fl6c|) . From f!7a|) . this implies /i~ 2 = 0, which is only 
possible for special values of s according to f!7bj) . It turns out that \v\ — 1 for the values 
of s that nullify 1/h. One concrete example has A = 1 and s = 7r/(2y / 3), while keeping an 
arbitrary nonzero B. Specifically, we have for this particular critical point of the differential 
system (Q: 



l/h(s) 



v(s) 



s=n/(2V3) 



s=tt/(2</3) 



B [1 + cos(2V3 s)] e 3 



=tt/(2V3) 



sm 



{V3s) 



t/(2V3) 



(8a) 



(8b) 
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The actual value for s at the critical point is nonphysical (because a is); what matters is that, 
for example, the Ricci scalar diverges there. For cosmic times just before the singularity, 
the functions v and 1/h can be expected to be slightly different from those given in (j7]). 
Moreover, the following corollary can be obtained from fl6cl) : 



,'\2 



3 1 -v 2 



(9a) 



- singularity 



where the suffix is interpreted as being arbitrarily close to the point with l/h(s) — dr/ds — 
corresponding to (IBa|) . The explicit solution ( ITal) can also be seen to satisfy (I9al) . In turn, 
JaJ gives for the vector- field energy density (I3"a|) the following result: 



?"vec 




I3\h 2 . 


singularity 



Pv 



3(E_ 



Planck 



) 2 H 2 



(9b) 



singularity 



with dimensional quantities in the middle expression. A final characteristic concerns the 
diverging vector-field equation-of-state parameter w vec = p ve c/r ve c and can be stated as 
follows: 



v/3(l -f 2 )/16 



-1 



(9c) 



singularity 



which, using (19bj) . can also be written with 3 A h 2 in the denominator. 

Further mathematical discussion of the final singularity flH]) is left to a future publication. 
Note that, strictly speaking, the qualification "final" is arbitrary, as the tensor- vector-scalar 
theory ([1]) is time-reversal-invariant and so is the differential system ([5]). 

For the present article, the relevant observation is that the numerical results of Fig. [1] can 
be interpreted as interpolating between the critical points fll]) and ([8]). Indeed, the particular 
combination 3 _1 h~ 2 (v) 2 / (1 — v 2 ) from the middle panel of Fig. [T]is seen to run between the 
values and 1, which are the corresponding values from (j2J) and fl9a|) . The numerical results 
for the ratio r vec / (3 A h 2 ) from the right panel of Fig. [1] show the same behavior, running 
between the values and 1 from fll]) and f!9bj) . respectively. Numerical results also match 
( l9cl) . but have not been shown explicitly in Fig. [TJ 

Including standard-matter, the final singularity is characterized by having p vec /A — > oo 
and pm/ A — > const > 0, in addition to having a diverging Ricci scalar R as mentioned 
before. This exotic behavior, just as that of Fig. [TJfor the pM = case, is, most likely, the 
result of the unusual properties of the vector-field energy density and pressure, which will 
be discussed in the next section. 



V. DISCUSSION 

Let us, first, elaborate on the remark of Sec. II III about the finite age of the type of model 
universe shown in Fig. [TJ For initial values of the standard-matter energy density that are 
not too large (compared to the value of A), one has an age of the order of 

^max - ^in ~ Qn -^Planck/ , (10a) 
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where the numerical coefficient o m depends on the vector- field boundary conditions at t = t{ D 
(cj n ~ 10 for the boundary conditions of Fig. [1]). From the analytic solution of the vector- 
field equation ({Sc]) with h replaced by l/v3, it is estimated that the dependence of c in on 
the initial values is only logarithmic, 



Ci In 



C 2 fin + C 3 Vh 



(E. 



Planck J 



1 V (t h 



dV (t) 



dt 



(10b) 
(10c) 
(lOd) 



t=t- m 



with positive constants c%, C2, C3 of order 1 and generic small values fllOcj) and (llOdj) . making 
for a positive and finite logarithm in ( UObj) . 

In a de Sitter spacetime with Hubble constant H^s, the Gibbons-Hawking tempera- 
ture 7gh_|23| effectively sets the scale of the initial vector-field perturbation by mode mix- 
ing Hii7so that \V (t in )\ ~ T GH ~ H dS ~ A^/Ep^ and |dVoM(tin)| ~ |#ds Vb(tto)|. 
The parametric dependence of fllOaj) is then given by 



, (GH-temperature) 
max in, perturbation 



(vector-field theory) 
dS 



In [(£pi anck ) 4 /A] £pi anck /VA 



(11) 



where t; n is the coordinate time of the low-frequency (long-wavelength) matter perturbation 
that breaks the original de Sitter symmetry and where the vector-field theory considered is 
the one given by (PQ). It is certainly possible that a result similar to f[TT]) can be obtained 
for other nonstandard matter fields, but it remains to determine precisely which types of 
matter fields suffice. 

For A — > + while keeping -Epi anc k fixed, the estimated lifetime ( ITT]) increases without 
bound, [t max — thjds — > +00. This behavior agrees with the naive expectation that the 
Minkowski solution of the A = theory remains effectively stable, even in the presence of 
the vector field V a (x). Note that the type of vector- field model considered with A > can 
still give an infinite-age solution (with Minkowski spacetime js| appearing asymptotically) 
if the energy-density function e(Qi) in (fTa]) is more complicated than the negative quadratic 
function ( lib)) ; see also later comments. 

The behavior found in Sees. II I II and [IV] differs from that of "normal" matter, which 



typically behaves according to the so-called cosmic-no-hair conjecture 23L l25l. l26|. Loosely 



speaking, the conjecture states that, with appropriate matter content, expanding universes 
that are not too irregular approach an eternal de Sitter universe (see also the recent pa- 
per 27| . which will be commented on below). For homogenous cosmological models, the 
cosmic-no-hair conjecture has been shown to hold 261 ) . provided the matter obeys both the 
strong energy condition (SEC) and the dominant energy condition (DEC). Recalling the 
succinct discussion of Ref. [28j] in terms of perfect fluids (here, for simplicity, specialized to 
the case of an isotropic pressure P), the SEC corresponds to having p + P>0Ap + 3P>0 
and the DEC to p > A P G [-p, +p\. 

The numerical solutions with only a standard-matter component (not shown here (l~4| ) 
agree with the expectations of the cosmic-no-hair conjecture. But not so for the solutions 
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with an additional nonstandard vector- field matter component: the model universe runs 
away from the de Sitter solution instead of towards it. The same conclusion holds for the 
rAi(r) = case presented in Fig. [TJ The numerical solutions display, in fact, a violation of 
the DEC on one count (P ^ [— p, +p] for p > 0) and a violation of the SEC on two counts 
(p + P < and p + 3P < 0). 

It is, however, clear that quite reasonable physical systems may display violations of the 
various energy conditions 29j, perhaps the least surprising being the violation of the SEC. 
In fact, SEC violation occurs already for a positive gravitating vacuum energy densit y (p y 



-Py > 0), which can result from underlying microscopic physical degrees of freedom 301.131 



The crucial question is whether the classical vector-field theory (Op) can be made into a 
consistent quantum theory. A related question is whether or not the vector-field theory 
([1]) can be shown to arise as an effective theory (see below for some remarks on infrared 
quantum effects [sill]). Obviously, the interest of the present article is only mathematical 
if the answer to both questions turns out to be negative. 

For completeness, it should be mentioned that the inapplicability of the cosmic-no-hair 
conjecture has also been discussed recently in the context of anisotropic inflationary models 
(cf. Ref. 23] and references therein). But the 'mild' behavior found in these anisotropic 
models 23] contrasts with the 'catastrophic' behavior resulting from the vector-field theory 
([TJ), where the isotropic model universe simply comes to an end (see Sec. [IV] and Refs. 15- 
|22j). Moreover, the theory considered in this article has a genuine positive cosmo logical 
constant A, not just a positive value of the scalar potential at a particular localized field 
configuration, which can make a difference for the nonlinear field equations and certainly 
does make a difference for the global spacetime structure 0]. 

In fact, the global structure of de Sitter spacetime has been argued to be responsible for 



a fundamental quantum instability through particle production [9l4l l| . The search is for a 



macroscopic description of the corresponding backreaction effects. Naively, our vector-field 
theory (CQ) appears to be ruled out, as the Hubble parameter increases due to the instability 
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m a way 



(details in the caption of Fig. [TJ see also the middle panel of Fig. 2 in Ref. 
reminiscent of what happens with an evaporating Schwarzschild black hole. However, the 
same type of vector-field theory can also give a decreasing Hubble parameter (see the top- 
right panel of Fig. 1 in Ref. [8J]), provided the quadratic energy- density function e(Qi) in 
(llap is replaced by a more complicated function El)]. Hence, if an effective vector- 

field theory is somehow relevant for the macroscopic description of backreaction effects from 
particle production js il], then the microscopic processes themselves select an appropriate 
macroscopic e-type function. 
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